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Abstract: We apply the QCD factorization approach to exclusive, radiative B meson 
decays in the region of small invariant photon mass. We calculate factorizable and non- 
factorizable corrections to leading order in the heavy quark mass expansion and next- 
to-leading order in the strong coupling constant. Phenomenological consequences for the 
, B if *7 decay rate and the B — > K*£^£^ forward-backward asymmetry are discussed. 



Radiative i^-meson decays provide an important tool to test the standard model of 
electroweak interactions and to constrain various models of new physics. The theoretical 
description of exclusive channels has to deal with hadronic uncertainties related to the 
binding of quarks in the initial and final states. For the decays B — > K*j and B — > 



D I K*£~^i , that we are focusing on here, this is usually phrased as the need to know the 

\ hadronic form factors for the B K''*^ transition, but there also exist "non-factorizable" 

strong interaction effects that do not correspond to form factors. They arise from the 
^ ■ matrix elements of purely hadronic operators in the weak effective Hamiltonian with a 

' photon radiated from one of the internal quarks. In Ref. Q we have computed these non- 

factorizable corrections and demonstrated that exclusive, radiative decays can be treated 
in a similarly systematic manner as their inclusive counterparts. As a result we obtain 
the branching fractions for B K*^ and B — > K*£~^i~ for small invariant mass of the 
lepton pair to next-to-leading logarithmic (NLL) order in renormalization-group improved 
perturbation theory. 

In the "naive" factorization approach, exclusive radiative B decays are described in 
terms of hadronic matrix elements of the electromagnetic penguin operator Oj and the 
semi-leptonic operators Og^io Q. These are parametrized in terms of the corresponding ten- 
sor, vector and axial-vector B K* transition form factors {Ti^2,3{q'^), Vig"^), ^o,i,2(^^))- 
Factorizable quark-loop contributions (Fig. Pa) with the four-quark operators Oi-q are 
taken into account by using "effective" Wilson-coefficients, Cj Cf, Cg Cf{q^), 
renormalized at the scale fi = mb- 
* Speaker. 

^Based on work together with M. Beneke and D. Seidel Q. 
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Figure 1: LO contributions to {'j*K*\Hcs\B). The circled cross marks the possible insertions of 
the virtual photon line. In (a) and (b) the spectator line is not shown. 



In order to include non-factorizable contributions as in Fig. ^ and Fig. |2|it is convenient 
to introduce generalized form factors Ti{q^) for the transition into a virtual photon B — > 
K*^* as follows, 



{Y{q,l,)K*ip',en\H,fi\B{p)) 



2T^{q^)e^^^^P-e:ppp'^-iT2{q^) [{Ml-ml,)e* 



{e*-q) ip^+p'n] 
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-{p^+p' 



K* 
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In the "naive" factorization approach these new functions reduce to %{q^) = Cf Ti{q^)+... 
Following the QCD factorization approach to exclusive B decays factorizable and non- 
factorizable radiative corrections are calculable in the heavy quark mass limit and for small 
photon virtualities (in practice < 4m^). 

At leading order (LO) in the strong coupling constant, the generalized form factors 

read 



T,{q')=CfT,{q')+Y{q') 



q 



2mb{MB + ruK* 



Viq'), 



%{q')=CfT2{q')+Y{q') 
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(2) 



The function Y{q'^), which is usually absorbed into C^{q^), arises from the quark loop in 
Fig. Pa. The last, "non-factorizable" term in T3{q'^) comes from the annihilation graph in 
Fig. ^ when the photon is emitted from the light spectator in the B meson (all other graphs 
are sub-leading in the l/mf, expansion). It introduces a new non-perturbative ingredient, 
namely one of the two light-cone distribution amplitudes of the B meson, 4>b,±{^), see 
||l|, ^ for details. Furthermore, for the considered values of g^, the recoil-energy of the 
out-going K* meson is large, and the seven independent B — > K* form factors can be 
described in terms of only two universal form factors Q, which we denote as S,±{q^) and 
^11 (g^) for transversely and longitudinally polarized K* mesons, respectively Q. 
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Figure 2: Non-factorizable NLO contributions to {'y*K*\Hcg\B). Diagrams that follow from (c) 
and (e) by symmetry are not shown. 



Figure 3: Factorizable NLO corrections to the B K* form factors. 



Factorizable next-to-leading order (NLO) form factor corrections are derived from 
Fig. ^ after the corresponding infra-red divergent pieces are absorbed into the soft universal 
form factors and ^y, see 0] for details. The non-factorizable vertex corrections (Fig. §c- 
e), are similar to the NLO calculation for the inclusive b — > S7* transition, and the result for 
the two- loop diagrams in Fig. |2|d+e are taken from Ref. Q. For the vertex corrections we 
chose a renormalization scale fj, = 0{mb). The non-factorizable hard-scattering corrections 
in Fig. ^+b and Fig. |^ involve the light-cone distribution amplitudes of both, B and K* 
mesons. (For = diagrams of this form have already been considered in Q, but using 
bound state model wave- functions, rather than light-cone distribution amplitudes.) Since 
in these class of diagrams the typical quark- and gluon-virtuality is of order y^AqcD^^ 
we chose a different renormalization scale fi' of that order. In principle, we also have to 
consider NLO order corrections to the annihilation graph in Fig. |^. However, since this 
term is suppressed by small Wilson coefficients C3 and C4 and numerically small already 
at LO, we have neglected these effects. Notice however, that the annihilation topology is 
numerically more important for B ^ decays ^. 

The B K*^ decay rate is proportional to the function |7i(0)p = |72(0)p. In order 
to study the effect of NLO corrections it is convenient to define a generalized exclusive 
"Wilson" coefficient C7 = 7i(0)/^_l(0). In Fig. ^ we have shown the ^-dependence of IC7P 
at leading order (LO), including only next-to-leading order vertex corrections (NLOi), 
and including all next-to- leading order corrections (NLO). As expected, the NLOi vertex 
corrections cancel the renormalization-scale dependence of the LO result to a great extent. 
(The hard-scattering corrections, arising at order reintroduce a mild scale-dependence.) 
Most importantly, we observe that the NLO corrections significantly increase the theoretical 
prediction for jC/p. Numerically, we have ICtInlo — 1-78 • ICyIlq- From this we predict the 
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branching ratio as 
Br(S 



K*-f) = (7.9 



+1.8\ 

i.eJ 



10 



-5 



1.6ps J V4.6GevJ ^ 0-35 ) 

Comparing with the current experimental averages [10| Br{B^ — > -^*^7)exp = (4.54 it 0.37) • 
10-^Br(5- ^^*-7)exp = (3.81±0.68)•10-^ and using the value C^(O) = 0.35 from QCD 
sum rules [0], we observe that the central value of the theoretical prediction overshoots 
the data by nearly a factor of two. (An equivalent analysis with similar conclusions can be 
found in Ref. §.) 

Possible explanations for this discrepancy 
are: i) new physics contributions (this is rather 
unlikely because of the good agreement between 
NLO theory and experiment for the inclusive 
counterpart, B Xs'j), ii) sizeable l/ruh power- 
corrections ( "chirally enhanced" corrections play 
a role for decays into light pseudoscalars [^] ; in 
our case, however, we expect a less dramatic 
effect), iii) an insufficient understanding of the 
B — > K* form factors (a fit to the experimen- 
tal data on the basis of our formalism yields a 
somewhat smaller value, '^±(0) = 0.24 it 0.06). 
A quantity that is less sensitive to the precise value of C±{q'^) is provided by the 
B K*i~^i~ forward-backward asymmetry ^fb- At LO the position of the asymmetry 
zero Qq is determined by the implicit relation 



0.25 



0.15 




0.05 



Figure 4: as a function of the renor 

malization scale /z, see text. 
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and does not depend on form factors at all |13]. As illustrated in Fig. |5| NLO cor- 
rections shift the position of the asymmetry zero from = 3.4^q5 GeV^ at LO to 
= 4.39;[JJ^ GeV^. (A shghtly different value = 3.94 GeV^ is found if one takes 
the complete form factors from QCD sum rules |11], instead of and the factorizable 
NLO corrections from Q). In any case, a measurement of the forward-backward asym- 
metry zero provide a clean test of the Wilson-coefficient Cg in the standard model with a 
rather small theoretical uncertainty of about 10%. 

In summary, we have shown that a systematic improvement of the theoretical descrip- 
tion of exclusive radiative B meson decays is possible. This is because in the heavy quark 
limit decay amplitudes factorize into perturbatively calculable hard-scattering kernels and 
universal soft form factors or light-cone distribution amplitudes, respectively. The next-to- 
leading order corrections increase the branching ratio for the decay B — > K*'y by almost a 
factor of two (which is at variance with the current experimental data if "standard" values 
for the soft form factors are used). They also shift the position of the forward-backward 
asymmetry in the decay B K*£~^i~ towards Qq = 4.2 it 0.6 GeV^ in the standard model. 
In this case the precision of the prediction is sufficient to test the Wilson coefficient Cg 
with only 10% theoretical uncertainty. 
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Figure 5: The FB asymmetry as a function of (left). The Wilson-coefficient Cg as a function 
of the FB asymmetry zero (right). The error band refers to a variation of all input parameters and 
changing the renormalization scale between mt,/2 and 2mf,. The dashed line is obtained from using 
the complete form factors from [jlll, see text. The grey band indicates the standard model value. 
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